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Strong gravitational lenses provide an important tool to measure masses in the distant Universe, 
thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch 
of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. 
However, the limiting factor in all of these studies has been the currently small samples of known 
gravitational lenses (~ 10^). The era of the SKA will transform our understanding of the Uni¬ 
verse with gravitational lensing, particularly at radio wavelengths where the number of known 
gravitational lenses will increase to ~ 10^. Here we discuss the technical requirements, expected 
outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We 
find that an all-sky (Stt sr) survey carried out with the SKA 1-MID array at an angular resolution 
of 0.25-0.5 arcsec and to a depth of 3 /iJy beam^* is required for studies of galaxy formation 
and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the 
SKAl is required for tests of dark matter and studies of supermassive black holes at high redshift 
to be made using gravitational lensing. 
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1. Introduction 

Gravitational lensing is the deflection of light from a distant background object {the source) 
by an intervening foreground mass distribution {the lens). If the surface mass density of the lens 
is sufficiently high, and if the source and lens are properly aligned, then multiple images of the 
background source are formed. The positions and distortions of these multiple images give impor¬ 
tant information about the mass distribution of distant lensing galaxies where velocity dispersions 
from spatially-resolved spectroscopy are either difficult or impossible. They also give a magnified 
view of the distant background source, probing a source population that is typically too faint to be 
detected with current instrument sensitivities. 

After the discovery of the first gravitational lens system B0957-I-561 (Walsh et al. 1979), in the 
Jodrell Bank 966-MHz quasar survey, the early history of strong lensing was dominated by VLA- 
based surveys, including the MIT surveys (Hewitt et al. 1992) of mainly steep-spectrum sources and 
the JVAS/CLASS surveys (Myers et al. 2003; Browne et al. 2003) of flat-spectrum radio sources. 
The majority of the currently known 36 radio-loud gravitational lens systems were found in this 
period. Further radio-based surveys became limited by the survey speed of high-resolution tele¬ 
scopes; because the typical separation of multiple images in gravitational lens systems, where the 
deflector is a massive galaxy, is ~ 1 arcsec, gravitational lens discovery programmes require large 
amounts of time on relatively oversubscribed interferometer arrays. Even with the recently up¬ 
graded JVLA and e-MERLIN, the survey speed is insufficient for large-scale searches. During the 
last 10 years, the advent of large-area optical imaging and spectroscopy surveys (e.g. SDSS) have 
found hundreds of gravitationally lensed star-forming galaxies (Bolton et al. 2008) and quasars 
(Inada et al. 2012). The advent of wide-field surveys in the sub-mm with Herschel (Negrello et al. 
2010) and the SPT (Vieira et al. 2013), coupled with the resolution and sensitivity of AEMA, is also 
expected to increase the numbers of strong gravitational lenses by a few hundred. In total, there are 
around 300-500 gravitational lenses known from galaxy, galaxy-group and galaxy-cluster lensing, 
with about 10 per cent having a radio-loud (-bright) background source. 

Gravitational lenses have been of extraordinary astrophysical value in many fields: for deter¬ 
mining the radial mass density profiles of galaxies (e.g. Koopmans et al. 2009); the direct detection 
of sub-galactic scale substructure predicted by dark matter models for galaxy formation (e.g. Veg- 
etti et al. 2012); measuring the mass of supermassive black holes in distant quiescent galaxies (e.g. 
Winn et al. 2004); measuring the stellar initial mass function in external galaxies by a combination 
of lensing and stellar population modelling (e.g. Auger et al. 2010); and by providing a one-step 
measurement of cosmological parameters, most notably the Hubble constant using time-delays 
(e.g. Suyu et al. 2010). Eurthermore, the lens also magnifies the background source by factors of 
~ 5-100, which gives unique information about the structure of high redshift objects. This allows 
a source population that is intrinsically faint to be observed and studied, for example, the faintest 
known radio sources are currently gravitationally lensed objects (e.g. Jackson 2011). Also, the 
magnification from the gravitational lens provides a high resolution (sub-kpc) view of galaxies at 
cosmologically distant epochs {z ~ 1-5; e.g. Deane et al. 2013). 

Most of these astrophysical applications have been limited by small-sample statistics, since 
only a minority of gravitational lenses are usually suitable for any particular study. However, we 
are on the verge of a huge increase in the number of gravitational lenses available for study in the 
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radio with the SKA and in the optieal/infrared with missions such as Euclid and the LSST. In this 
review, we show that a lensing survey with the SKA, coupled with the synergy provided by Euclid 
and the LSST, will open up a new parameter space, both in terms of statistics and in the discovery of 
rare and valuable gravitational lenses for studying galaxy formation, the high redshift Universe and 
cosmology. This review builds upon the Radio All-sky SKA Lensing (RASKAL) survey proposed 
by Koopmans et al. (2004) during the first SKA Science Workshop. 

2. Discovery space of SKA: The gravitational lensing statistics of continuum sources 

Surveys such as those envisaged with the SKA will be capable of increasing the sample of 
radio-loud gravitational lens systems from a few tens to tens of thousands relatively straightfor¬ 
wardly. In order to determine the actual numbers that are expected to be found, we need to evaluate 
the lensing probability, T(Zi), of a given radio source, 

f cd t 

t{zs)= n{z)Ox—dz, (2.1) 

Jo dz 

where Zs is the source redshift, n{z) is the number density of the lenses. Ox is the cross-section 
of the gravitational lenses, and cdt is the path length to the background source (e.g. Turner et al. 
1984). By integrating the lensing probability over all of the sources in the parent sample, we can 
determine the total number of sources that are gravitationally lensed. If we know the number counts 
of the parent sample as a function of flux-density and redshift, then we can calculate the number of 
gravitational lenses expected within a flux-limited survey. To zeroth order, gravitational lensing is 
a rare event; a high-redshift object has a ~ lO^^-lO^^ chance of lying within the Einstein radius 
of an intervening galaxy-scale mass distribution (depending on the source population redshift and 
magnification bias). For cluster lensing the probability is even smaller, at about 10 Therefore, 
surveys of a few hundred million high redshift objects should in principle give access to a few hun¬ 
dred thousand lenses, although recognition difficulties will mean that the actual number discovered 
will be a factor 10-100 less than this (see Section 3). 

In order to estimate the lensing potential of the SKA surveys, we first require a model that 
describes the lens and source populations. To calculate the number density of lenses, we have 
used the light-cones generated by Kitzbichler & White (2007) coupled with semi-analytic pre¬ 
scriptions (De Lucia & Blaizot 2007) from the Millennium simulation (Springel et al. 2005). We 
have assumed that each potential lens has a velocity dispersion that is related to its luminosity via 
the Faber-Jackson relation, with a normalisation that gives for an L* galaxy a velocity dispersion 
of a* = 190 kms^^ We have then scaled the L* luminosity using the lensing statistics from the 
CFASS survey (Myers et al. 2003); this is done by requiring that an initial sample of radio sources 
with flux densities > 1 mJy, and differential source counts with a slope of —2 (McKean et al. 2007) 
produces 22 gravitational lenses, as observed by CFASS (Browne et al. 2003). This has been done 
by assuming the CFASS selection criteria, namely that selected gravitational lenses should have 
flux ratios < 10:1 and an image separation > 300 mas. To first order, this calibration removes the 
effects of cosmic variance. The SKADS database (Wilman et al. 2008) was used to provide the 
population of background radio sources. This population is a best-guess extrapolation of known 
radio surveys down to the /xJy level (e.g. Condon et al. 2012), and will almost certainly be su¬ 
perseded once the SKA surveys begin in earnest. The exact nature of the very faint (< 0.1 /xJy) 
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population is important, because the magnification provided by lensing (factors 5-100) yields many 
more gravitational lens systems in populations with a steep source count, a phenomenon used to 
advantage in siih-mm surveys (e.g. Negrello et al. 2010). 

In total, 4.56 x 10^ radio sources from a 1 square degree area of SKADS and with a 1.4-GHz 
flux density > 1 /iJy were extracted and randomly placed in the Millennium area. Any galaxy 
within 1 arcsec was considered to be a potential lens; the lensing status was then evaluated for each 
such pair. The results from the simulation are shown in Figure 1. In the 1 square degree of SKADS 
used, 440 radio sources with a total radio flux density > 1 /iJy were found to be gravitationally 
lensed and have image separations > 0.3 arcsec. Of these lensed sources, 49 have total flux- 
densities > 10 /iJy, and 10 have total flux-densities > 50 jxiy. These calculations clearly show 
that for a fixed observing time it is most efficient to carry out a shallow, wide-field survey for 
gravitational lenses with the SKA. As we will demonstrate in Section 3, a survey with an rms 
of ~ 3 /xJy beam^' (2 year duration), and a detection threshold of about 15a for lensing, yields 
about ~5 lenses per pointing (0.5 square degrees) using the full 250 dish SKAl-MID array. These 
systems will contain lensed objects typical of the radio population at a slightly fainter flux density 
level of a few tens of /iJy, in other words, approximately evenly distributed between lensed AGN 
and lensed starbursts. 

The simulations also show the gravitational lenses that may be available from surveys with the 
SKA2, which will be up to a factor of 5-10 deeper than with SKAl-MID. These are potentially far 
more exotic objects, with the lensing galaxies at redshifts 1^, extremely faint, very high redshift 
background sources, and in very large numbers despite the unfavourably flat faint radio source 
count of SKADS; the lensing rate approaches ~ 1:150 due to the high redshifts (median z ~ 3) 
of the background source population. The details of these samples should be taken with caution 
because they rely on uncertain extrapolations. In particular, the likelihood of a large population 
of radio sources at z > 10 is low. What these simulations do illustrate, however, is the unique 
power of a gravitational lensing survey to probe any faint, high redshift radio source populations 
that may exist, if we know the detailed properties of the lensing galaxies. The latter is currently 
unknown, but is an area of very active study that will progress greatly over the next 5-10 years. For 
example, a deep (sub-jUJy) small area continuum survey with the SKAl-MID array would directly 
detect this population of magnified and exotic high redshift objects, accessing the type of objects 
that will be seen routinely with the SKA2. Alternatively, if the predicted number of lenses is not 
found, then this would mean that the faint radio source population has a much different luminosity 
function than is currently expected. However, this assumes that the challenges related to the lens 
identification and confirmation process are well understood, which is typically the case for lens 
surveys are radio wavelengths where the lens selection function (resolution, sensitivity, flux-ratios; 
e.g. Myers et al. 2003) is well defined. 

The major uncertainties in the number of gravitational lenses that will be potentially available 
for study in any SKA survey are the extrapolation of the radio source counts to very faint flux den¬ 
sity levels and high redshifts, and the detection efficiency. The former uncertainty is negligible at 
the level of a few tens of /iJy (which is most relevant for an SKAl survey), but is very probably a 
factor of a few at the 1-ftJy level, and likely an order of magnitude for the z > 5 radio source pop¬ 
ulation. This uncertainty dominates the estimates for any SKA2 survey. The detection efficiency 
has an uncertainty that is likely to be of a factor of a few, although more detailed blind simulations 
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Figure 1; The results from the SKA gravitational lensing simulations. Upper left. Lens galaxy redshift 
versus source galaxy redshift. Upper middle: Lens galaxy redshift versus source galaxy flux density. Upper 
right: Lens galaxy magnitude versus source galaxy flux density. Lower left: Lens galaxy Einstein radius 
versus source galaxy flux density. Lower middle: The predicted differential number counts normalised to 
Euclidian from the SKADS simulation. Lower right: The total integrated number counts as a function 
of flux density for all sources (dark blue), the lensed sources (red), the lensed sources using the unlensed 
flux-density (green) and the lensing incidence (cyan). 

are needed to reduce this uncertainty. 

3. Gravitational lens candidate identification and confirmation techniques 

The main observational challenge will be the recognition of genuine gravitational lenses among 
the large numbers of false positives corresponding to radio sources with intrinsic structure; the key 
criterion is angular resolution, although frequency coverage and image sensitivity also play a role. 
Known galaxy-scale gravitational lenses have image separations between ~ 0.3-4.5 arcsec, with 
the typical image separation being ~ 1 arcsec (e.g. McKean et al. 2005). From our SKA lens¬ 
ing simulations, the predicted image separation distribution for sources with a total flux-density 
> 50 /f Jy has a range between 0.3 and 1.4 arcsec. Therefore, an angular resolution of between 
0.25-0.5 arcsec is required for identifying most of the lens candidates for follow-up imaging. It is 
for this reason that gravitational lensing surveys will only be possible with the SKAl-MID array. 

To simulate the expected imaging quality for a gravitational lens search, we have taken the 
array configuration for SKAl-MID and generated mock gravitational lensing data based on an 
extended star-forming galaxy and a compact AGN; for both types of object, we have used total 
source flux-densities of 1, 0.25 and 0.05 mJy. From Figure 2, we see that sensitive observations 
(Briggs weighting; Robust = 0) with SKAl-MID have sufficient sensitivity (~ 3 /f Jy beam^^) and 
angular resolution (~ 0.5 arcsec) to detect galaxy-scale gravitational lenses with extended sources. 


5 

















Strong gravitational tensing 


]. P. McKean 


We find that it is possible to achieve an angular resolution of 0.25 arcsec (for Robust = 2), but as 
expected, this increases the image rms by a factor of ~ 5, leading to the fainter extended radio 
sources having a surface brightness below the detection limit. Nevertheless, observations at Band 
2 (0.95-1.76 GHz) with SKAl-MID can be used to select gravitational lensing candidates. An 
all-sky (in sr) survey with SKAl-MID, reaching a sensitivity of 3 /xJy beam^^ could potentially 
find ~ 3 X 10^ gravitational lens systems with image separations > 0.3 arcsec. 

The resulting recognition problem (candidate to confirmed lens) will be considerable, but can 
be eased by adding additional information. The first is the large bandwidth to be used in the 
SKAl-MID survey, which will mean that some spectral index information will be available to 
discriminate between multiple images of the same object and, for example, a non-lensed source 
consisting of a flat-spectrum core and steep-spectrum jet or lobe emission. In addition, follow-up 
observations of candidates at Band 5 (4.6-13.8 GHz) will give a more robust measurement of the 
specttal index of the components and provide higher angular resolution imaging of the candidates 
(0.06-0.20 arcsec; Robust = 0 weighting) to determine if the surface brightness of the images is 
conserved, as is required by gravitational lensing. The second piece of information will come from 
surveys at other wavelengths, in particular those of Euclid and the LSST. Cross-correlating the SKA 
gravitational lens candidates with the high resolution optical data will determine whether there is a 
massive lensing galaxy at the expected position (see Figure 3) and provide a photometric redshift 
estimate of the gravitational lens. Even if neither the SKAl-MID and EuclidESST surveys are 
capable of detecting gravitational lenses on their own, experience from CLASS shows that the use 
of two surveys at different wavelengths (radio and optical) gives a very much greater discrimination 
between multiply-imaged sources and those that display intrinsic structure. 

Using the predicted number counts and image separation distribution from our SKADS sim¬ 
ulation, we have estimated the number of gravitational lenses that should be straightforwardly 
identified from the SKAl-MID survey presented above with minimum follow-up. Here, we only 
consider those lensed images with separations > 0.5 arcsec and total flux-densities > 1 mJy. We 
find a conservative estimate of about ~ lO'* new radio-loud gravitational lens systems would be 
found (~ 1 per 3 square degrees), which is a factor of over 10^ more than is currently known. We 
note that an initial pilot survey, that covers the same sky area as the Euclid deep fields (15000 deg^) 
and to an rms of 6 /iJy beam^', would conservatively identify around 2000 strong gravitational 
lenses in 3 months observing time. An SKAl-MID array that is only 50 per cent complete would 
still be useful for finding gravitational lenses, providing the sensitivity and length of the baselines 
remains unchanged and a wide area survey is carried out. The situation will improve dramatically 
with the SKA2, where the improved sensitivity, particularly on the long baselines will allow the 
potentially large population of ~ 10^ faint lensed radio sources to be straightforwardly identified 
in the survey imaging data. 

4. Searches for gravitationally lensed spectral line sources 

An alternative method to efficiently find new gravitational lens systems with the SKA is to 
select candidates via magnification bias; here a specific property of the source results in very steep 
number counts. For example, Negrello et al. (2010) successfully used this technique to select 
a sample of strong gravitational lenses by using a 500/im flux density limit and exploiting the 
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Figure 2: Simulations of an SKA 1-MID snapshot observation of an extended radio galaxy (top) and an 
AGN (bottom) with total flux-densities of 1, 0.25 and 0.05 mJy (left middle, right middle, right) using 
Briggs weighing Robust = 0 (0.5 arcsec beam size) and for a total flux-density of 1 mJy (left) with Robust 
= 2 (0.25 arcsec beam size). 



Figure 3: High resolution HST optical imaging of the gravitational lens CLASS B0631H-519 with the MER¬ 
LIN 1.7 GHz radio contours shown overlaid. The optical data only detect the two lensing galaxies, whereas 
the radio data are only sensitive to the emission from the background source. Together, these two datasets 
clearly show that CLASS B0631H-519 is a gravitational lens (York et al. 2005). 


steep slope of siih-mm number counts to obtain a ~ 100 per cent reliable strong gravitational lens 
sample. Gonzalez-Nuevo et al. (2012) extended this technique to the high-luminosity end of the 
sub-mm luminosity function. The HI mass function is expected to broadly follow its local form of a 
Schechter function (e.g. Abdalla et al. 2010). Therefore, a wide-held SKA HI survey, as proposed 
by Abdalla et al. (2010), should also be amenable to magnification bias at the high-mass end of the 
HI mass function. 

For the detailed formalism used for calculating differential magnification probability distribu- 


7 
























Strong gravitational tensing 


]. P. McKean 


tions, /?(/!,z)d/i, for a strong gravitational lens with magnification ft > 2, see Perrotta et al. (2002, 
2003). It is expected that the high magnification tail has the form f>(ft,z) = for some func¬ 

tion a{z), regardless of the lens population. The a{z) function depends on the nature and evolving 
number density of the gravitational lenses. Blain (1996) and Perrotta et al. (2002, 2003) consider 
several options for galaxy lens populations. In the latter, the mass spectrum follows the Sheth & 
Tormen (1999) formalism. As for the continuum case (Section 2), we assume a singular isothermal 
sphere density profile for fhe lens populafion and assume a non-evolving model, buf normalised fo 
fhe Peroffa el al. predicfions al an arbilrary redshifl of z = 1. 

A more significanf uncerlainly in fhe lensing predicfions is fhe maximum magnificafion caused 
by fhe finile source sizes. We follow Perrofla ef al. (2002) in spanning fhe range of plausi¬ 
ble maximum magnificafions wifh ft < 10 and ft < 30, implying source characteristic radii of 
~ 1 — 10 kpc; note that if the assumed source sizes are larger than this then the maximum 
magnifications will be lower, but it is expected that galaxies will be systematically more compact 
at higher redshifts (e.g. Gunn & Gott 1972) so our assumed sizes are likely fair. Several higher 
magnification events are known (e.g. Swinbank et al. 2010), but we conservatively neglect this 
more extreme population. Note that the surface density of strong gravitational lenses increases 
with maximum magnification. The choices of maximum magnification here is more conservative 
than that of, for example, Blain (1996) who used ft < 40, and Lima et al. (2010) who used ft < 100. 

Figure 4 shows the results of applying this strong gravitational lensing formalism to a non¬ 
evolving HI mass function (Zwaan et al. 2003), and to the evolving mass function of Abdalla 
et al. (2010). Selecting gravitational lens candidates as those with HI masses of > 12 M* in an HI 
spectral line survey, should yield a sample of between 0.5-5 gravitational lenses per square degree, 
but with a ~ 100 per cent selection efficiency and the key advantage that the source redshifts are 
known. Further details, and additional calculations for the Dark Energy Grism survey with Euclid, 
are presented by Serjeant (2014). 

5. Science capabilities with SKA 

In this Section, we review the main astrophysical results that will come from the strong grav¬ 
itational lensing programme of the SKA. These will be a combination of results that will exploit 
the large statistical sample of gravitational lenses, or will use smaller samples of rare gravitational 
lenses for investigating galaxy formation and cosmology. 

5.1 Galaxy formation studies 

The hierarchical model for galaxy formation predicts that structure in the Universe formed 
through the mergers of smaller mass systems, resulting in the tapestry galaxies, clusters of galaxies 
and large-scale structure that we observe today. A crucial result of this model is that the stellar and 
gas components of these structures should be embedded within an extended dai'k matter halo (e.g. 
Navarro et al. 1996). The form of the radial mass-density profile of the resulting dark matter haloes 
and the level of low mass substructure within them are both sensitive to the galaxy formation model 
and the energy of the dark matter particle. 

Over the last decade, gravitational lenses have been used to place important constraints on 
the mass distributions for galaxies at cosmological distances and have tested these predictions 
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Figure 4: The HI mass function (black line) at redshifts z = 0.5,1,2, assuming no evolution from the z = 0 
mass function (top) and assuming an evolving mass function (Bottom). The no-evolution lens population 
is shown in blue, and the singular isothermal sphere lenses of Perrotta et al. (2002, 2003) are in pink. The 
hatched regions show the range of maximum magnification jU < 10 and jj. < 30. Note that at > 12 M* the 
observed population is dominated almost entirely by gravitationally lensed HI galaxies. 


from hierarchical galaxy formation models. At radio wavelengths, high resolution imaging of 
gravitationally lensed extended radio sources at mas-scales with VLBI have constrained the inner 
(within 5-20 kpc) radial mass density profiles for a handful of lens galaxies (e.g. Cohn et al. 2001; 
Wucknitz et al. 2004), finding fhaf fhey are close fo isofhermal (i.e. 7 = 2 where p(r) oc r^l'). 
However, fhe mosf significanl resulfs have come al opfical wavelengfhs, where high resolution 
imaging wifh fhe HST of fhe Einslein rings and arcs of extended blue sfar-forming galaxies, coupled 
wifh stellar kinemafics (velocify dispersions and velocify fields) for fhe lens, have defermined radial 
mass densily profiles for a homogenous sample of 58 early-fype galaxies at z ~ 0.2 (Koopmans et al. 
2009). Again, these results find fhaf massive galaxies have on average isofhermal densily profiles 
(7 = 2.09 ±0.02). 

An imporfanl conslrainl on galaxy formalion models from gravilafional lensing comes fhrough 
fhe defection of low mass subslrucfure around more massive galaxy-scale haloes. The cold dark 
mailer model for galaxy formation predicls a subslrucfure mass function, dn oc (e.g. Die- 

mand ef al. 2007), and a mass fracfion in subslrucfure of < 1%. However, fhe lack of observed 
safellifes found around our own Milky Way suggesls fhaf fhe mass funclion may be much flaller 
lhan Ihoughl, which could be due lo our Milky Way being a unique case, fhe subslruclures being 
Ihere in fhe expected abundance, bul are dark, or due lo a differing physics of fhe dark mailer (e.g. 
warm dark mailer; self-inleracling dark mailer). Gravilafional lenses can be used lo lesl galaxy 
formalion models on Ihese small-scales because fhe surface brighlness disfribulion of fhe observed 
images depends on fhe form of fhe polenlial and fhey can be (very) sensitive lo local perturbations 
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in the mass model due to small elumps of dark matter in the parent halo (or along the line-of-sight). 
There are two main applieations of this teehnique. The first is gravitational imaging (Koopmans 
2005; Vegetti & Koopmans 2009), whieh uses the information eontained within extended ares to 
deteet the presenee of individual substruetures and measure their mass properties. This teehnique 
has been used to deteet substruetures with masses as low as 2 x 10^ M© out to redshift z = 0.881 
from infrared adaptive optics imaging at ~ 65 mas resolution (Vegetti et al. 2012). The second 
method flux-ratio anomalies (e.g. Metcalf & Madau 2001; Dalai & Kochanek 2002), which 
occur when the flux-density of one or more point images are different from what is predicted by a 
smooth mass model. This method is sensitive to the total amount of substructure at the Einstein ra¬ 
dius of the lens, and can be used to statistically detect a substructure population down to 10“^ Mq. 
Here, observations at radio wavelengths are important because the flux-ratios are less effected by 
microlensing from stars in the lensing galaxy (Koopmans et al. 2003) and extinction due to dust is 
not an issue. 

The SKA, when combined with optical telescopes like Euclid and the LSST, will revolutionise 
studies of the mass distributions of galaxies by increasing the overall sample sizes and will allow 
measurements to be made for samples over a wide range of galaxy mass, type, redshift and envi¬ 
ronment. The main results that will come from the SKA strong gravitational lensing sample are as 
follows. 

• The main constraints on the radial mass-density profiles of sfrucfure will come from high 
resolution imaging faken in Band 5 (4.6-13.8 GHz) wifh fhe SKAl-MID array, which can 
pofenlially be coupled wifh HI velocify field measuremenfs of lenses al low redshifl (z < 0.4) 
using specfral line observations in Band 2; combining fhe consfrainfs from fhe lensing and 
dynamics breaks fhe degeneracies befween mass models fhaf are inherenf in bofh mefhods 
individually (e.g. Barnabe el al. 2009). The complelion of fhe SKA2 will allow a much 
larger population of lenses al higher redshifl lo be sludied in Ibis way, which will determine 
whelher Ihere is any evolulion in fhe mass properlies of galaxies wifh epoch and provide a 
key observational conslrainf lo galaxy formalion simulations. 

• The very precise measurement of the projected mass within the Einstein radii of galaxies 
(~ 10 kpc) from strong lensing can be used, with optical and infrared measurements, to 
independently determine the stellar mass and place constraints on the stellar initial mass 
function (e.g. Auger et al. 2010). This requires both the redshift of the lens and source to 
be known, and the modelling of the SKAl-MID array data to determine the Einstein radii 
of the lenses. The SKA strong gravitational lenses will provide a large statistical sample to 
investigate variations of the stellar IME with galaxy type and mass when coupled with the 
optical/IR imaging provided by Euclid and the ESST. 

• The SKAl-MID array has the potential to provide important constraints at the low end of 
the dark matter halo mass function using strong gravitational lensing, which would indepen¬ 
dently test models for the dark matter particle. The main requirement for the gravitational 
imaging method is the angular resolution of the data (source structure and signal-to-noise 
ratio are also important). Eollow-up observations with SKAl-MID in Band 5 of those lenses 
with extended radio sources will give an angular resolution of 0.03-0.1 arcsec, similar to the 
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resolution provided by the HST and ground based adaptive optics telecopes, and so should 
give detections of substructures down to ~ 10^ Mq for a large sample of lens galaxies. How¬ 
ever, the most important constraints to the dark matter halo mass function will come from 
observations on VLBI-scales since the mas angular resolution will be sensitive to structures 
at the ~ 10^ Mq level (see Figure 5), a mass regime where various models for the dark 
matter particle strongly differ. Therefore, an SKAl-MID array that is capable of VLBI ob¬ 
servations with, for example, the EVN, is required for this unique science goal of the strong 
gravitational lensing programme to be achieved. 

• It will also be possible to place constraints on the substructure mass function with the SKAl- 
MID array using those gravitational lenses with compact (point) emission from AGN (see 
Figure 2). It is expected that there will be ~ 300 four image gravitational lens systems from 
the SKAl-MID survey that can be used for this; a factor of about 50 more than are currently 
known. This will require the precise measurement of the flux-ratios of the point-like images; 
this will be done through monitoring the systems with the SKAl-MID Band 5 to achieve 
the desired angular resolution and to determine if the fluxes have been changed by intrinsic 
variability and the gravitational lensing time-delay (the latter can be used to determine the 
cosmological parameters, see below). 

• The sensitivity of the SKA will allow the routine detection of highly de-magnified core lensed 

images that are predicted to form at < 100 pc from the centre of the lensing galaxy (for shal¬ 
lower than isothermal lensing potentials). Such images probe the central mass distribution 
and can be used to constrain the central density profile where dark matter is insignificant, but 
the influence of both the central supermassive black hole (SMBH) and the central stellar cusp 
are important (e.g. Mao et al. 2001). Understanding the central regions of galaxies is impor¬ 
tant because they are directly relevant to the question of how the SMBH influences (and is 
influenced by) the process of galaxy formation. This central question has been debated since 
the discovery of the relation between the SMBH mass and the bulge stellar velocity disper¬ 
sion on much larger scales (Ferrarese & Memtt 2000). However, the detection of core-lensed 
images is extremely rare due to their very high de-mganification ( 10 ^-10 ^) relative to 

the other lensed images. There is currently only one confirmed case (Winn et al. 2004). 

Targeted observations with SKAl-MID should allow the detection of core-lensed images 
in about 100 systems, compared to the 20-30 accessible to e-MERLIN, provided that an 
angular resolution of about < 100 mas can be achieved, as is expected from observations 
at Band 5 with the SKAl-MID. The large frequency coverage of Band 5 is also required 
to confirm that the emission is from the central image and not the lensing galaxy (see Eig- 
ure 6). With the SKA2, measurements of the central-image properties of thousands of lens 
systems will be possible, and hence will provide a census of central potentials in many dif¬ 
ferent lensing galaxies as a function of redshift, mass and environment. With these sensitiv¬ 
ities, non-detections become just as interesting as detections because they typically indicate 
SMBH masses greater than the M-o relation, unusually steep central density profiles, or 
both. Moreover, in combination with VEBI studies, as could be done if the SKA operated 
as part of a VEBI array, it may be possible to detect splitting of the central image on mas- 
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Figure 5: Simulations of the effect low mass dark matter haloes have on arc surface brightness distributions 
for a 10® Mq (left), 10^ M© (middle) and 10^ M© (right) substructure. The top row of images show the 
full gravitational arc, detected on VLBI scales with an SKAl-MID+VLBI array, and the bottom row show 
a zoomed in view of the region most affected by the substructure. The 10^ Mq and 10^ Mq substructures 
can be seen directly in the data, whereas the 10® M© substructure would need to be detected using the 
gravitational imaging technique. 

scales due to the combined effect of the black hole and central stellar cusp; in these cases all 
degeneracies can be broken and both components’ masses can be measured separately. 


5.2 Studies of the lensed source population 

Although the SKA will transform our understanding of the high redshift Universe, the sen¬ 
sitivity and angular resolution of the telescope can be increased by orders of magnitude (at no 
additional cost) through using the magnification provided by gravitational lenses. In this respect, 
almost any high-redshift object science application with the SKA could benefit from the use of nat¬ 
ural telescopes, providing a model for the lens (the optics) is well understood; there are now several 
sophisticated methods available for the reconstruction of the distorted sources behind gravitational 
lenses (e.g. Wucknitz 2004; Vegetti & Koopmans 2009; Hezaveh et al. 2013). 

Here we briefly summarise some of the principle applications of the SKA strong gravitational 
lensing sample for studying the background source population. 

• The increase in sensitivity (factors of 5-100) will allow a population of objects that could not 
otherwise be observed to be detected and studied in detail. For example, the faintest known 
continuum radio source (S 5 ghz ~ 1 MJy; Jackson 2011) is gravitationally lensed (see Figure 
7). The first immediate result will be to study the parent population of sources that will 
be detected with the SKA2 using the SKAl-MID array. The form of the parent population 
luminosity function is sensitive to the number of gravitational lenses that will be detected 
in the deep continuum surveys to be caiTied out with SKAl-MID. For example, about 440 
gravitational lenses should be detectable in the SKAl-MID ultra deep survey (1 deg^; rms 
0.05 /f Jy; see Prandoni & Seymour 2014). Furthermore, the magnification provided by the 
gravitational lenses (~ 10 ) will enable radio sources with star-formation rates of > 1 Mq yr^^ 
at z ~ 3-4 and > 5 Mq yr^^ at z ~ 6 to be detected. These sources would directly probe the 
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Figure 6: The gravitational lens system J1632-0033 is the first and only gravitational lens system to have 
a central lensed images for a galaxy-scale system. Images A and B are of the redshift 3.42 lensed radio 
sources, whereas image C is the core lensed image. The spectra of the three images are identical are high 
frequencies, since propagation effects (free-free absorption within the lensing galaxy) are less severe (Winn 
et al. 2004). 


radio luminosity function in a regime that is around an order magnitude lower in luminosity 
than would be possible without the magnification from gravitational lensing. 

• The magnification will be particularly important for studying spectral line sources for which 
the large observational bandwidths will not increase the detectability, but will increase the 
volume that can be searched. An example is presented in Figure 7, where the most dis¬ 
tant radio spectral line in emission is shown (water maser at z = 2.64; Impellizzeri et al. 
2008). Water maser systems are expected to be highly abundant at high redshift (McKean 
et al. 2011) and could be used to measure the masses of supermassive black holes in distant 
galaxies through observations of the maser line kinematics over time (e.g. see Kuo et al. 
2013 for studies at low redshift). Coupled with high angular resolution imaging (see be¬ 
low) of highly magnified water maser lines, it could also be possible to measure geometric 
distances for radio sources at moderate redshifts (z ~ 0.7-1) and test models for dark en¬ 
ergy. In addition, water maser emission is known to be highly variable and so could also 
provide an independent measurement of dark energy through measuring the time-delay due 
to gravitational lensing (see below). High resolution observations would be possible using 
the SKAl-MID array Band 5 (probing z ~ 0.6-4) for detection and variability monitoring, 
and with an SKAl-MID array that is part of a global VLBI array for imaging. Also, it may 
be possible to detect the HI absorption forrest from the Epoch of Reionization by targeting 
the population of high redshift objects that have had their intrinsic flux density boosted by 


13 






Strong gravitational tensing 


]. P. McKean 





Figure 7: (left) JVLA imaging of the highly magnified radio source (images A, B, C, D) behind the lensing 
cluster SDSS J1004+4112 (Jackson 2011). (middle) The spectral line emission detected from the most 
distant known water maser galaxy, MG J0414+0534 at redshift 2.64 (Impellizzeri et al. 2008). (right) The 
CO (2-1) Einstein ring from the z = 4.12 quasar PSS J2322+1944 at 0.3 arcsec resolution with the VLA 
(right) and the reconstructed source (left) after lens modelling, which shows a rotation disk of gas (Riechers 
et al. 2008). 


gravitational lensing. Such observations would use the SKAl-LOW array. 

• The magnification provided by the lens increases the solid angle of the background source, 
which allows the object to be observed at high angular resolution. Such observations will be 
most useful for studying the morphology of radio sources at high redshifts and the structure 
of their various thermal and non-thermal components on pc to sub-kpc scales; e.g. star- 
formation, atomic and molecular gas distributions and AGN activity. An example of what 
can potentially be done is also shown in Figure 7, where the CO molecular gas distribution 
of a redshift ~ 4 quasar (Riechers et al. 2008) has been determined after correcting for the 
effect of the lens. This will require a high frequency capability (< 30 GHz) to investigate 
high redshift galaxies (z > 3), and the full collecting area of the SKA2 to study the formation 
of galaxies at the Epoch of Reionization (z ~ 6-10). 


5.3 Studies of cosmology 

Gravitational lenses can be used to measure distances at cosmological redshifts, and are there¬ 
fore potentially important for cosmology. The most direct application of gravitational lensing for 
cosmology is the measurement of the Hubble constant, which can be achieved by monitoring gravi¬ 
tational lenses with background sources that are variable in flux-density. As the images result from 
multiple different light-paths, the measurement of time delays between image variations gives a 
measurement of the difference in angular diameter distance, and if the redshifts of the source and 
lens are known, this in turn measures the Hubble constant, together with a secondary dependence 
on other cosmological parameters such as and Ha. The first Hq measurements were conducted 
on gravitational lenses discovered in the radio, mainly using the VLA in the 1990s (e.g. Biggs 
et al. 1999; Fassnacht et al. 2002). Subsequent Hq determinations were predominantly made us¬ 
ing gravitational lens systems containing radio-quiet, optically-selected quasars, mainly due to the 
availability of dedicated small optical telescopes, for example, recent measurements of time-delays 
have been dominated by the COSMOGRAIL collaboration (e.g. Eulaers et al. 2013). 
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The major systematic problem with gravitational lens determinations of Hq is that, as well 
as the time delays, accurate mass models of the lensing galaxies and the surrounding masses are 
needed. This process is subject to complex degeneracies (e.g. Schneider & Sluse 2013), and needs 
careful modelling and treatment of the systematics (e.g. Suyu et al. 2010) to achieve robust results. 
Nevertheless, two well-modelled gravitational lenses in which the systematics are under control 
give errors on Hq (dominated by systematics) of about 5% {Hq = kms^^ Mpc^^; Suyu 

et al. 2013), and further progress is likely in the next few years. Since Hq is not determined by 
CMB measurements unless assumptions are made about other cosmological parameters (such as 
zero curvature or w = —1), very accurate standalone measurements of Hq are capable of improving 
other cosmological constraints; for example, the dark energy figure of merit can be improved by 
factors of several (Linder 2011), and the more reliable gravitational lensing results already provide 
a constraint of w = — 1.521 q2q at the 68 per cent confidence level (Suyu ef al. 2014). 

The second major cosmological applicafion of gravifafional lensing uses fhe sfafisfics of fhe 
numbers of gravifafional lenses or fhe properfies of gravifafional lenses in well-selecfed samples. 
These are sensitive bofh fo fhe numbers and masses of fhe pofenfial lensing galaxies al parlicular 
redshifls (and hence fo galaxy evolufion), fo fhe disfribufion and luminosify funcfion of fhe sources, 
and fo fhe cosmological world model. Thus, we can use any samples in which fhe lenses and 
sources are well undersfood fo consfrain cosmology, or use a world model, fogefher wifh lens 
surveys fo investigate galaxy evolution. In fhe early days, gravifafional lensing sfafisfics were used 
fo confirm fhe exislence of A-like terms in fhe world model (e.g. = 0.691 q27 ; Chae el al. 

2002 ); laler, gravifafional lensing sfafisfics were used fo provide evidence for (lack of) evolution in 
elliptical galaxies (Chae & Mao 2003). Recenlly Oguri el al. (2012) have argued lhal lens sfafisfics 
consfrain more robuslly, = O.VOj^QQy (slal.) (sysl.), and are consislenl wifh a w = —1 
equation of slale, allhough fhe laller requires fhe combinafion wifh olher cosmological probes. 

The main cosmological applications of fhe SKA gravifafional lens sample are as follows. 

• Our simulations indicate lhal fhe SKA will provide a much more powerful probe of gravi- 
lafional lenses in fhe high-redshifl universe (see Figure 1). On fhe mosl optimistic assump¬ 
tion, lhal of isolhermal polenlials in low-mass, high-redshifl lensing galaxies, vasl numbers 
of high-redshifl lenses will be produced. The likelihood lhal Ihis assumplion is false, and 
lhal such galaxies will be more dark-malfer dominafed and hence less efficienl gravifafional 
lenses, ilself implies lhal sfafisfics of an SKA lensing survey will conlain information aboul 
mass dislribulions in high-redshifl galaxies. Improvemenfs in fhe cosmological world model, 
in parlicular in w, are likely fo resull from fhe discovery and followup of a large number of 
new lens syslems. In parlicular, fhe lensing slalislics of samples from CLASS and fhe SDSS 
are based on samples of jusl 13 and 19 gravifafional lenses, respectively. Increasing such 
samples by orders of magnilude, will resull in fhe uncerfainlies being dominated by fhe sys- 
lemalics related to fhe assumed galaxy evolufion model. 

• Increasing fhe samples of gravifafional lenses by a faclor of > 10^ will allow us fo monilor 
gravifafional lenses for lime delays and be more seleclive in using Ihose gravifafional lenses 
for which fhe syslemalics can be well conlrolled. This will require fhe AGN-dominaled 
systems fo be observed and monitored wifh fhe SKAl-MID array in Band 5 to provide belter 
dala for gravifafional lens modelling and fo delermine lime-delay lighl curves (observalions 
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at 8.46 GHz at 0.2 arcsec resolution have been suceessfully applied in the past). Constraints 
on w are then possible by combining the gravitational lensing time-delay data (which are 
mainly sensitive to Hq) and the CMB data from Planck (e.g. Suyu et al. 2013). Time-delay 
measurements will become more powerful both with numbers, and with a larger range of 
redshifts of gravitational lenses for which information can be extracted, as the statistical 
uncertainties are dominated by the handful of currently known systems that are appropriate 
for such an analysis. 

• Increasing samples by factors of > 10^ will also yield rare and cosmologically important 
objects. The ideal lens for cosmology would be one in which two sources are present at 
different redshift (e.g. Collett & Auger 2014), but in which one is variable, such as a radio 
quasar. This allows an angular distance measurement without problems of mass degeneracy, 
providing kinematical information for the gravitational lens is known; a few such lenses give 
w essentially free of systematic errors (see Schneider 2014 for discussion of the systemat- 
ics). Such systems will be a factor 500-1000 rarer than “standard” radio lenses, because the 
lensing probability on any one line-of-sight is 0.1-0.2%. This is within the likely yield of the 
SKA2 surveys. 

6. Summary 

The SKA era will transform our understanding of the Universe and strong gravitational lenses 
will play an important role in investigating galaxy formation and evolution over cosmic time, the 
properties of dark matter and the cosmological world model. This review has demonstrated that 
the SKAl-MID array is capable of detecting ~ 10^ new radio-loud gravitational lenses, providing 
there is sufficient sensitivity on 0.25-0.5 arcsec baselines and a frequency coverage between 1-14 
GHz. A more conservative estimate suggests that ~ lO'* should be straightforwardly identified in 
an SKAl-MID array survey al ~ 3 /ilJy beam^' sensilivily and wilh a 371 sr sky area, wilh Ihe main 
uncerlainlies being dominaled by Ihe luminosily funclion of Ihe background source populalion. 
Coupling Ihe SKAl-MID dala wilh opfical and infrared surveys, and furlher high resolulion radio 
imaging, will be needed lo idenlify Ihe remaining gravilalional lenses in Ihe expecled wide-area sky 
surveys. An inilial pilol survey, mafched lo Ihe sky area covered wilh Euclid should find several 
Ihousand new gravilalional lenses from only 3 monlhs of observing wilh Ihe SKAl-MID array. 
Such a pilol survey would lesl Ihe lensing slalislics and lens idenlificalion melhods. Finally, we 
note lhal il is imporlanl for several of Ihe unique science goals of Ihe SKA slrong gravilalional 
lensing programme lhal Ihere is Ihe VLBI capabilily of Ihe SKA. 
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